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von  Karman  constant  in  the  model  were  examined  with  results  shown. 
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ABSTRACT 


Specification  of  parameters  defining  the  mean,  as  well  as  the 
turbulent  properties ,  of  the  over  water  surface  layer  is  examined .  A 
barocllnic  boundary  layer  model  is  considered  using  Cardone's  modifi¬ 
cations  for  the  marine  environment. 

Data  fields  of  pressure  and  temperature  were  selected  as  input, 
and  original  ship  observations  were  utilized  to  verify  the  predicted 
surface  parameters .  The  effects  of  stability  changes ,  isobaric  curva¬ 
ture,  thermal  wind,  and  the  value  of  the  von  Karman  constant  in  the 
model  were  examined  with  results  shown. 

Iteration  criteria  used  in  the  model  were  evaluated  with  respect 
to  operational  application.  Winds  calculated  by  the  model  at  the  height 
of  19.5  meters  were  incorporated  into  the  current  FNWC  wave  height 
program,  and  the  results  compared  with  observed  wave  data.  Turbulent 
parameters  describing  the  properties  of  the  index  of  refraction  (C^)  were 
also  computed  and  evaluated  with  respect  to  laser  propagation  criteria. 
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I.  INTRODUCTION 


Parameterization  of  the  boundary  layer  over  the  sea  surface  has 
been  receiving  increased  attention  in  recent  years.  Fortunately  much 
of  the  theory  developed  for  low  level  turbulence  over  a  land  surface  is 
also  applicable  over  the  oceans .  But  the  air-water  interface  does  pre¬ 
sent  some  unique  problems.  For  instance  a  fluid  surface,  rather  than  a 
fixed  one,  results  in  mobile  surface  roughness  elements.  Because  of 
the  presence  of  a  fluid  surface ,  any  parameterization  requires  a  knowl¬ 
edge  of  the  processes  which  transfer  heat  and  momentum  energy  across 
the  boundary  and  ultimately  influence  such  variables  as  wave  generation 
and  index  of  refraction. 

The  purpose  of  this  study  was  to  examine  and  evaluate  several 
features  of  a  recently  formulated  marine  boundary  layer  model  which  was 
designed  for  operational  use.  The  model  is  a  two- layer  baroclinic  model 
for  the  marine  boundary  layer  and  the  formulation  examined  was  that 
developed  by  Cardone  (1969).  Cardone  based  his  formulation  on  an 
earlier  two-layer  model  due  to  Blackadar  (1965a).  Features  of  interest 
in  the  model  are  related  to  the  physics  of  the  planetary  boundary  layer 
and  include  the  effect  of  atmospheric  stability,  baroclinlcity ,  coriolis 
acceleration  and  changes  in  the  surface  roughness  parameter  on  the 
pressure  driven  surface  wind. 


13 


The  scope  of  the  study  was  to  examine  the  model  in  the  context  of 
its  operational  usage.  As  such,  it  was  a  case  study  in  which  the  data 
represented  pressure  and  temperature  observations  from  sequential 
12-hour  periods  when  a  large  cyclonic  pressure  system  dominated  the 
Central  North  Pacific  Ocean.  This  particular  synoptic  situation  provided 
the  necessary  relations  between  the  scalar  fields  to  examine  features 
associated  with  barocllnlclty ,  stability  and  wind  speed  categories. 

Because  these  data  represent  scalar  fields  which  had  to  be  de¬ 
fined  by  some  analysis  scheme  before  being  used  as  input  to  the  two- 
layer  model,  the  differences  which  arise  due  to  procedures  in  the  external 
field  specifications  are  considered .  The  effect  of  isobaric  curvature  on 
the  surface  wind  fields  is  examined  and  this  effect  is  compared  to  the 
changes  associated  with  the  physics  in  the  model.  The  differences 
which  may  arise  solely  because  the  fields  are  defined  by  an  objective 
versus  a  subjective  analysis  are  considered . 

Finally,  expressions  describing  the  boundary  layer  are,  by  their 
very  nature,  empirical.  As  such,  they  contain  empirical  constants  and 
the  von  Karman  constant  (k)  is  one  whose  value  has  recently  been  ques¬ 
tioned.  The  importance  of  a  suggested  change  from  k=0.40  to  0.35  was 
considered. 
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n.  ASPECTS  OF  A  PLANETARY  BOUNDARY  LAYER 


A.  SURFACE  BOUNDARY  LAYER 

1.  Basic  Structure  of  the  Boundary  Laver 

Recent  advances  In  the  study  of  atmospheric  turbulence 
have  led  to  a  greater  confidence  In  describing  the -distribution  of  wind 
and  temperature  in  the  atmospheric  boundary  layer  over  a  solid  surface. 
This  progress  has  arisen  largely  through  the  application  of  the  Monin- 
Obukhov  similarity  theory  which  has  been  supported  by  recent  observa¬ 
tional  studies.  Cardone's  contribution  was  primarily  to  extend 

Blackadar's  two- layer  neutral  baroclinic  boundary  layer  model  over  fixed 

«• 

terrain,  and  then  to  apply  this  model  to  the  boundary  layer  over  a  fluid 
surface.  The  major  difference  oetween  these  boundaries  is  that  the  sea 
surface  changes  its  character  as  the  wind  blows  over  it. 

In  general,  Blackadar's  model  treats  the  boundary  layer  as 
being  formed  by  both  turbulent  friction  and  the  corlolis  acceleration.  All 
processes  are  considered  to  be  three  dimensional,  and  the  layer  is  char¬ 
acterized  by  a  density  stratification.  It  is  conveniently  separated  into 
the  three  regions  depicted  in  Figure  1 .  In  this  delineation  the  important 
consideration  is  the  specification  of  which  is  the  turbulent  transfer 
coefficient  for  momentum. 
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Fig.  1.  Diagram  of  the  planetary  boundary  layer 


a.  Viscous  Sublayer  (less  than  1  cm) 

In  this  near  surface  layer  the  stress  Is  supported  by 
viscosity  (aerodynamlcally  smooth  flow),  and  the  mean  velocity 
increases  linearly  according  to 


(1) 


where  v  Is  the  kinematic  viscosity  and  U#  is  the  friction  velocity. 

b.  Surface  Layer  (up  to  50  meters) 

This  is  the  constant  stress  layer  in  which  the  coriolis 

parameter  can  be  neglected  and  K  increases  linearly  with  height.  For 

m 

neutral  stratification  the  logarithmic  law  is  valid 
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where  k-0.4  Is  the  von  KArmln  constant. 

c.  Ekman  Layer  (50-600  meters) 

This  layer  Is  characterized  by  the  stress  decreasing 

with  height  and  K  constant  throughout  the  layer, 
m 

In  formulating  the  boundary  layer  model,  Cardone  (1969) 
considered  the  validity  of  the  "log  law"  for  the  surface  layer  over  water, 
which  has  been  experimentally  documented  down  to  near  surface  limits 
(Phillips,  1966;  Paulson,  1967).  He  also  assumed  that  the  process  of 
wave  generation  did  not  seriously  affect  the  surface  layer  since  under 
neutral  and  steady  conditions  a  logarithmic  shear  zone  should  be  estab¬ 
lished.  Accepting  the  log  profile  as  being  valid,  Cardone  was  able  to 
apply  much  of  tho  surface  layer  turbulence  theory  developed  over  land 
to  over-water  conditions .  However  Cardone  did  take  Into  account  the 
action  of  the  waves  when  specifying  the  effective  roughness  parameter 
for  the  sea  surface. 

2.  Similarity  Theory  for  the  Surface  Layer 

The  Monln-Obukhov  "similarity  theory"  has  led  to  consider¬ 
able  progress  In  the  specification  of  the  distribution  of  wind  and  tem¬ 
perature  in  the  boundary  layer.  The  basis  of  this  theory  is  that  near 
the  ground  there  exist  velocity  (U#),  length  (L),  and  temperature 
scales  which  are  essentially  Invariant  with  height.  When  the  principle 
varlebles  such  as  temperature  (T),  wind  (U),  and  height  (z)  are  expressed 


at  nondimentional  fractions  of  theta  quantities ,  a  series  of  nondimen- 
sional  equations  can  be  formed  which  are  of  general  validity  in  the  sur¬ 
face  boundary  layer  (Lumley  and  Panofsky,  1964).  Recent  analyses  of 
the  oceanic  wind  data  (e.g.  Phelps  and  Pond,  1971)  have  indicated  that 
the  Monln-Obukhov  theory  is  valid  for  over-water  application.  These 
results  support  Cardona's  inclusion  of  the  Monin-Obukhov  stability 
criteria  in  the  boundary  layer  model. 

Similarity  theory  predicts  that  a  universal  relation  should 
exist  for  the  nondlmensional  wind  shear 


dU 

U*  dz 


(3) 


and  the  nondlmensional  temperature  gradient 


0*  dz 


ft 


(4) 


where  0t  and  0U  are  functions  of  the  dimensionless  ratio  z/L.  Here 
9  is  the  potential  temperature,  9*  is  the  scaling  temperature  defined  by 


H 


9*  kU*  Cp  0a  (5) 

and  L  is  the  Lettau- Monln-Obukhov  stability  length  defined  by 


-u2  Cd  Pa  T 
kgH 


(6) 


The  heat  flux  (H)  is  not  easily  measured,  so  a  modified  stability  length 
was  defined  by  Panofsky  (1963) 


Kh  U#T  du/dz 

L'  '  *hL  '  ^TL  '  k» 


(7) 
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The  assumption  that  L’  is  Independent  of  height  Implies  that  a^  is  con¬ 
stant,  which  In  turn  implies  similarity  between  wind  and  temperature 
profiles . 

If  L'  is  a  valid  scaling  length,  similarity  theory  predicts 


that 


U*  dz 


J*u  (z/L*) 


_l_  I® 
0*  dz 


0U  (z/V) 

ah 


(8) 

(9) 


Integration  of  these  formulae  (Panofsky,  1963)  yields 


Uz  -  Hi.  j^ln  z/zq  -  V(zA*)j 

^  ♦  ““ “  «/*o-  V(z/L')j  , 


where  the  relation  between  V  and  0U  is 


V(z/L‘) 


■/ 


l/L'  l-0„  <{> 


i 


d  C  ,  (  -  z/L' . 


(10) 


(11) 


(12) 


Paulson  (1970),  using  the  KEYPS  formula,  integrated  (12)  by  parts  to 
obtain  the  following  expression  for  V 


V  -  l-0u-3lnffu  +  21n  )  +  2  tan”  Vu  -  “5“  +  ln 


2  9  ’  ~  P'y  2  •  '  2 


Within  the  scope  of  an  operational  application,  these 
expressions  for  the  surface  layer  would  be  used  to  determine  the  wind 
speed  given  U*  ,  or  to  determine  U*  given  the  wind  speed  at  a  particular 
level.  For  this  reason  the  partial  derivatives  in  (7)  are  replaced  by 
expressions  involving  height  differences  and  air-sea  temperature 
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differences.  This  is  achieved  by  substituting  the  following  expression 
,  obtained  by  combining  (9)  and  (11) 


for 


dz 


ae  ^u  4?)  (®a "  ®s ) 

"  t  [ln(z»/*o)->t<(zVL')]  * 


(14) 


where  za  Is  the  height  at  which  0a  Is  measured  (chosen  as  10  meters  by 
Cardone) . 


Substituting  (14)  and  (8)  Into  (7)  yields 

U?  8  [in  (If)-  4»(ZaA')] 
k2g  (8a  -  0S) 

Equation  (10)  Is  then  solved  for  U* 


(IS) 


U* 


k  U 


m 


In  (2m/z0)  -  V(*m/L') 


(16) 


where  zm  Is  the  height  at  which  the  wind  speed  is  measured  (chosen  as 
19. 5  meters  by  Cardone).  Equations  (15)  and  (16)  can  then  be  solved 
for  U*  given  the  wind  speed  at  some  level  (Uz)  and  the  air-sea  tempera¬ 
ture  difference . 

3.  Consideration  of  the  Drag  Coefficient  (Cz)  and  Roughness 
Parameter  (zo) 

The  surface  roughness  parameter  and  the  Influence  of  thermal 
stratification  have  received  considerable  attention  In  the  definition  of 
the  relation  between  the  nondimenslonal  drag  coefficient  (Cz)  and  the 
wind  speed  at  a  certain  height  (usually  10  meters).  The  need  for  such  a 
relation  arises  from  the  "bulk  aerodynamic"  formula  for  estimating 
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momentum  flux.  This  equation,  deduced  from  dimensional  analysis,  has 
the  form 


-  P_  uj  -  p.  c.  u*  , 


(17) 


where  T  is  the  surface  stress.  From  this  expression  the  drag  coefficient 
le  defined  as 


(lb) 


2 

Rearranging  and  squaring  (10)  yields  the  following  expression  for  the 
drag  coefficient  at  a  level  z  in  terms  of  zQ  and  stability 


n2 


In  z/zn-  4»(7L') 


(19) 


which  was  the  form  considered  by  Cardone. 

Observations  have  indicated  that  over  the  sea,  zQ  is  depend¬ 
ent  on  U*.  This  dependence,  however,  changes  as  the  wind  speed 
changes.  For  the  case  of  light  winds  (aerodynamically  smooth  flow), 
the  drag  coefficient  Cz  decreases  with  increasing  Wind  speed  (Kraus, 
1966)  and  zQ  has  the  form 


V 

2o  =  1.09  U*  * 

For  aerodynamically  rough  flow,  C_  increases  with  increasing 
speed  according  to 

a  U  l 

Zo  g 

where  (a)  is  a  constants  0.035  (Charnock,  1955).  A  summary 


(20) 

wind 


(21) 
of  field 
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observations  Indicates  that  the  change  between  rough  and  smooth  flow 
occurs  at  wind  speeds  of  5-6  meters/second. 

Cardone  proposed  the  following  simple  expression  for 
covering  a  wide  range  of  wind  speeds  (both  rough  and  smooth  flow) 

^1  2 

Z0  =  “■  +  c2  u;  +  c3  .  (22) 

The  constants  were  chosen  so  that  Cz  was  a  minimum  for  a  value  of 
6  meters/second  at  the  10  meter  height  under  neutral  conditions ,  and 
above  this  speed  the  relation  approximated  Charnock's  results.  For  U* 
in  meters/second  and  zQ  in  meters,  the  resulting  expression  was 

zQ=  6.84  x  10"5/U*  +  4.28  x  10-3U*  -  4.43 x  10“.4  (23) 

Other  representations  of  C  versus  U  have  been  published . 

z  z 

A  few  of  these  appear  in  Figure  2  along  with  Cardone' s  approximation. 

It  is  clear  from  the  scatter  of  curves  in  Figure  2  that  the  specification 
‘of  the  drag  coefficient  was  a  primary  decision  in  Cardone' s  formulation. 
This  choice  of  C_  versus  U„  will  be  examined  during  the  discussion  on 
results  from  this  investigation. 

With  the  equations  for  stability  length  (15),  friction  velocity 
(16),  and  surface  roughness  parameter  (23),  the  surface  layer  wind  dis¬ 
tribution  is  completely  described  within  the  framework  of  the  similarity 
theory.  These  three  equations  can  be  solved  simultaneously  for  U* 
using  iteration  techniques  with  the  external  parameters  described  below 


as  input. 


2  4  6  8  10  12  14 
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Pig,  2.  Various  suggested  forms  of  the  drag 

coefficient  (Cz)  versus  wlndspeed  (U  ) 
at  10  meters 


B.  PLANETARY  BOUNDARY  LAYER  j  I 

.  \ 

l.l  Background  \ 

Ih  developing  a  model  for  the  entire  boundary  layer, 


con¬ 


sideration  has  to  be  given  to  the  pos  sible  external  parameters .  In 
Cardone's  formulation  these  parameters  were  the  geostrophipwind  (G), 
the  air-sea  temperature  difference  (0a  —  ) ,  the  coriolis  parameter  (f), 

and  the  thermal  wind  (r)  ^  He  considered  previously  suggested  empirical 
relations  between  the  geos  trophic  drag  coefficient  (U*/G)  and  the  sur¬ 


face  Rossby  number  (Lettau,  1959) 


R  =  7s- 

°  ^z0 


(24) 


where  G  is  the  surface  gdostrophic  wind'  speed.  Empirical  relations 

i  * 

have  also  been  suggested  between  the  cross  isobar  flow  angle  and 


V 


\ 


Blackadar  (1965a)  developed  a  two- layer  boundary  layer 


lyer 

I  I  ' 

model,  consisting  of  a  Prandtl  surface  layer  and  an  Ekman  layer,  which 

I 

was  verified  on  the  basis  of  the  above  relations. 

\ 

Cardone's  formulation  was  closely  related  to  Blackadar' s 
work.  In  particular)  the  height  of  the  surface  layer  (h)  was  specified 

explicitly  in  terms  of  the  external  parameters 

'  I 

Bo  G 


h  = 


(25) 


where  Bn  wai  a  nondimensional  constant  equal  to  3. Ox  10  ^  as  deter- 

1 

j 

mined  by  Blackadar. 


By  applying  proper  boundary  conditions  and  requiring  that 
continuity  of  wind  and  wind  shear  exist  at  level  h ,  the  following  solu- 

i 

I 

tions  for  U*/G  and  4>0  (along  with  the  solution  for  zQ)  completely 
define  the  surface  boundary  layer  wind  distribution  under  neutral 


conditions . 


*■[ 


2  k  Bq  sin 


*«] 


k^T  sin  (n/4-4»0) 


to  Bo  Ro 


Implicit  in  such  an  approach  is  that  the  value  of  Km  at  the  internal 

I 

boundary  between  the  surface  layer  and  the  Ekman  layer  is  given  by 


K  =| 
m 


kU*  G 


2.  Cardone's  Extension  of  Blackadar's  Two- Layer  Model 

i 

The  spccess  of  Blackadar's  two-layer  model  in  predicting 

the  essential  characteristics  of  the  neutral  surface  boundary  layer,  along 

\  l 

with  its  relative  ease  of  application,  suggested  its  extension  to  the  non¬ 
neutral  marine  boundary  layer.  Cardone  extended  the  diabatlc  model  to 

1  I 

over- water  flow  by  regarding  z0  as  a  function  of  U*  . 

The  eddy  viscosity  in  the  constant  stress  layer  for  the  non¬ 
neutral  case  was  expressed  by  Cardone  in  terms  of  the  nondimens ional 


wind  shear  as 


Kin  = 


kU*z 

0U(Z/L') 
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Since  Kjjj  Is  constant  with  height  above  h,  its  value  in  the  Ekman  layer 
is  then 


kU*  B0G 

,  Km  “  {0U  (Z/L‘)  * 

The  wind  speed  at  h  was  given  by 


(30) 


u*  r  . 

uh  -  —  m  b0r0  -  «*»  (h/L')  . 


(31) 


Internal  boundary  conditions  applied  at  level  h  appear  in  Figure  3  and 

dW  3  n 

follow  from  the  fact  that  *r —  always  makes  an  angle  of  — —  to 

OZ  4 

the  ageostrophic  wind  (W)  and  that  is  parallel  to  U  in  the 

a  z 

constant  stress  layer. 


Fig.  3.  Wind  vector  diagram  applicable  at  level  h. 

By  requiring  continuity  of  wind  shear  at  h,  as  well  as  directional 
continuity,  Cardone  was  able  to  reduce  the  problem  (analogous  to  (26) 
and  (27))  to 


u* 

G 

[2  k  B0  sin2  ipo  0U  (h/L')]  / 

(32) 

u*  _ 

G 

sin 

In  B0R0  - 

( n/4  -  V0) 

(h/L') 

(33) 
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u* 

The«e  equations  are  solved  simultaneously  for  —  and  q»  from  the 

G  o 

external  parameters  G,  f,  and  (9a-fls).  In  addition,  a  dimensionless 


stability  parameter  was  defined  as 

L* 


B0G 


f  L' 


(34) 


which  represented  the  ratio  of  the  surface  layer  h  to  the  stability 
length  L'  (Blackadar  and  Ching,  1965). 

The  relationships  determined  by  Cardone  for  U*/G  and 
as  functions  of  RQ  and  L*  appear  In  Figures  4  and  5.  For  all  surface 
Rossby  numbers,  U*/G  increases  and  decreases  with  decreasing 

stability.  In  general,  stable  conditions  have  a  more  pronounced  effect 
on  the  surface  layer  wind  characteristics  than  unstable  conditions .  In 
comparing  the  results  from  the  model  with  observations,  Cardone  noted 
that  U*/G  increased  by  roughly  30%  in  unstable  conditions,  with  no 
marked  tendency  to  decrease  at  great  instability,  and  decreased  roughly 
70%  in  stable  conditions.  These  features  in  the  observed  results  will 
be  discussed  later. 

The  effects  of  baroclinicity  were  also  included  by  assuming 
that  the  geostrophic  wind  was  a  linear  function  of  height  (Blackadar, 
1965b).  Blackadar  (1965a)  introduced  a  parameter,  p,  which  was  defined 
at  level  h  by 

p  -  -  TjJT/ V)  ["  B0R0  -Wth/I.  )]  .  (35) 


The  requirement  of  continuity  of  wind  sheer  for  the  barocllnlc  case  took 


the  form 

+  •  <3w 

Figure  6  Illustrates  the  geometrical  Implication  of  (36).  From  the  Ekman 
layer  solution  of  the  ageostrophlc  wind,  the  following  equation  can  be 
obtained 

p  Hr|h  -  V?  v  iwhi  . 


u. 


4 


dz 


where  p*  Is  a  dimensionless  parameter  defined  by 


P* 


BoG 

*u  <  Vl-  ) 


J  <„<"/«••> 

T  2kU*BcG 


['nBoRo- 


(38) 


Fig.  6.  Geometrical  relationships  resulting  from  the 
requirement  of  the  continuity  of  wind  shear  at 
the  internal  boundary  level  h  (Blackadar,  1965b). 

Next  the  thermal  wind  vector  was  replaced  by  its  dimensionless 


I  i 


d  z 


magnitude 


r 


f 


(39) 


d  G 

where  the  direction  of  relative  to  the  eurface  geos  trophic  wind 

wae  implied  by  angle  n  In  Figure  6.  Still  another  dimensionless 
quantity  was  defined 

r'  -  rB  f  In  B  R  -  UJ  (h/L’)]  (40) 

o  L  0  0 

and  (3S)  was  then  rewritten  as 

ph»f'|  ‘  r  G>  m 

By  applying  the  law  of  cosines  to  each  of  the  two  smaller  triangles  in 
Figure  6  that  have  M  for  one  side,  the  following  relations  were  determined 


M2  - 

K|2.2-  gV 

(42) 

2 

s  * 

(1  ♦  2p’  +  2p'2  ) 

(43) 

2 

q  - 

(1  +  r'2  -  2r*  cos  H  ). 

(44) 

0  is  the  angle  between  the  geostrophlc  wind  and  the  thermal  wind  (see 


Figure  6).  Then  the  law  of  sines  was  applied  to  yield 

sin  4  ■  P'/s  (45) 

sin  y  ■  sin  H  (46) 

d 

The  following  relations  apply  to  the  large  triangle  (ABC)  in  Figure  6 . 

uh  ■  |Wh  |  sin  (  •  ♦  y  )/sin  UIQ  (47) 

U2  -  G2  ♦  |  wh  |2  -  2  G  |wh  |  cos  (  •  ♦  y  )  (48) 

|wh|2  -  U2  ♦  G2  -  2  UhG  cos  mo  .  (49) 

From  equations  (31),  (45),  (46),  (47),  (48),  and  (49),  the  following 


equations  describing  the  surface  boundary  layer  characteristics  were 
obtained 
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Tan  III 


sin  (  *  +  >  ) 
(s/q)-  co*  (  •  ♦  t  ) 


U  -  *  G  q  »ln  (  •  ♦  y  ) 

*  •  tin  III  flnBR  -  MICVL')]  * 


(50) 

(51) 


O  ‘  0  0 

By  Including  (23)  and  (15),  Cardona  obtained  a  aystem  of  equations  that 

could  be  solved  for  U*  and  Ul  from  the  input  parameters  G,  f, 

o 


(6a-9g),  j^~|  ,  and  fl  •  These  equations  were  calculated  from  a 

knowledge  of  sea  level  pressure,  surface  air  temperature,  sea  surface 
temperature,  and  latitude. 

3.  Some  Characteristics  of  the  Marine  Surface  Boundary  Layer 
A  primary  effort  of  Cardone's  work  was  to  determine  how 
well  the  planetary  boundary  layer  model,  using  only  large-scale  synoptic 
parameters  as  Input,  specified  the  surface  boundary  layer  wind  distribu¬ 
tion  over  the  sea  surface.  Cardonc  considered  several  studies  relating 
the  geostrophlc  wind  ratio  Uz/G  and  the  air-sea  temperature  difference. 
The  results  of  an  extensive  study  by  Carstensen  (1967)  at  Fleet  Numerical 
Weather  Central  (FNWC)  appear  in  Figure  7  along  with  earlier  results 
due  to  Bleeker.  The  ratio  of  surfaco  wind  (V)  to  geostrophlc  wind  (G) 
is  tabulated  for  each  degree  of  the  reported  air-sea  temperature  differ¬ 
ence,  and  the  median  value  for  each  Interval  is  plotted  as  an  x  In 
Figure  7.  The  validity  of  the  results  was  limited  to  the  air-sea  tempera¬ 
ture  difference  range  -4  to  ♦PC.  In  order  to  Include  observed  winds 
in  the  surface  pressure  analysis,  FNWC  uses  a  constant  value  for  V/G 
of  0.78  and  a  turning  angle  of  15°  toward  lower  pressure.  This  selection 


could  produce  significant  error. 
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It  -  I.  "0 

Fig.  7.  Reported  wlnd-geostrophlc  *lnd  ratio  as  a 
function  of  air-sea  teaperature  difference 
(Carstensen,  1967) 
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Pig.  8.  Two  lajer  nodal  predictions  of  tho  dependence 
of  the  ratio  of  the  19*5  neter  wind  speed  to 
the  geostrophlo  wind  speed  upon  air-sea 
teaperature  difference  and  thermal  advectlon 
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Cardona's  model  can  be  used  to  compute  a  wind  speed  at 
19. 5  meters.  Figure  8  Illustrates  the  variation  of  the  ratio  of  the  19. 5 
meter  wind  and  surface  geostrophlc  wind  with  the  air-sea  temperature 
difference.  The  ratio  varied  less  In  unstable  conditions  than  In  stable 
conditions . 

The  effects  of  barocllnlclty  also  appear  In  Figure  8  and 
these  effects  lncreaso  as  stability  decreases,  with  the  ratio.  Ujg  /G , 
Increasing  with  cold  advectlon  and  decreasing  with  warm  advectlon. 
Mendenhall  (1967)  also  studied  the  Influence  of  barocllnlclty  and  sta¬ 
bility  on  the  veering  of  the  wind  in  the  planetary  layer.  He  concluded 
that  over  mid-latitude  oceans  the  diurnal  variation  of  frictional  veering 
of  the  wind  with  height  was  negligible,  but  that  lapse  rate  and  baro- 
cllnlcity  were  important.  Results  obtained  by  Cardone,  using  the  plan¬ 
etary  model,  showed  qualitative  agreement  with  Mendenhall's  conclu¬ 
sions.  Cardone's  results  also  suggested  that  barocllnlclty  was  more 
significant  than  stability  In  determining  the  veering  of  the  wind  In  the 
planetary  boundary  layer.  For  both  cases  in  situations  of  strong,  cold 
advectlon  the  actual  wind  may  back  with  height,  especially  in  conjunc¬ 
tion  with  unstable  stratification. 

Additional  results  on  the  influence  of  baroclinicity  came 
from  a  study  by  Clarke  (1970),  who  found  a  tendency  for  the  boundary 
layer  wind  to  back  Instead  of  veer  with  height  in  unstable  conditions . 
He  stated  that  this  could  have  been  due  to  barocllnlclty  or  to  accelera¬ 
tion.  He  noted  that  It  could  have  occurred  if  was  systematically 

o  z 
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negative  in  convective  conditions  up  to  a  level  of  about  0.2  "j — 
Instead  of  being  positive  as  in  the  normal  case.  This  latter  interpreta¬ 
tion  was  based  on  the  observation  that  the  wind  appeared  to  be  sub- 
geostrophlc  in  the  convective  layer  and  reached  geostrophic  values  only 
after  a  rapid  increase  through  the  overlying  inversion. 
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HI.  COMPUTATIONS  AND  NUMERICAL  PROCEDURES 


A.  EXTERNAL  WIND  FIELD  SPECIFICATIONS 

All  computations  were  made  with  data  defined  on  a  subset  of  the 
FNWC  63X63  grid  which  covered  the  North  Pacific  Ocean  (Figure  9). 

Ehta fields  included  sea  level  pressure,  sea  surface  temperature,  and 
surface  air  temperature  at  each  of  315  grid  points  spaced  381  km  apart 
at  60*N.  The  necessary  map  factor  was  included  for  a  1:30  million 
polar  stereographic  projection. 

The  corlolis  term  appears  as  an  external  parameter  and  therefore 
poses  a  limitation  on  the  lowest  geographical  latitude  where  the  geo- 
strophic  relation  would  be  valid.  It  has  already  been  established  that 
the  dynamic  boundary  layer  was  formed  by  the  surface  stress  and  the 
corlolis  acceleration.  Pushistov  (1970)  concluded  that  manifestations 
of  the  nonlinearity  in  the  dynamics  of  the  boundary  layer  and  the  appear¬ 
ance  of  a  perceptible  ageos trophic  component  of  velocity  in  the  free 
atmosphere  were  to  be  expected  near  latitudes  of  7-10°.  Sheppard 
( 1970)  suggested  20°N  as  a  limiting  latitude  for  application  to  boundary 
layer  models.  A  limiting  value  of  15°  latitude  was  used  in  this  program 
as  the  limit  on  geostrophic  flow.  For  latitudes  below  15°,  the  corlolis 
parameter  would  be  kept  constant. 

1.  Geostrophic  Wind  Computations 

The  geostrophic  wind  components  were  computed  from  the 
surface  pressure  field  (P)  as  follows 
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a 


V 


_£P 

by 

_J_P 

dx 


(52) 


(53) 


using  a  centered  explicit  finite  difference  technique.  The  density  was 
determined  from  the  equation  of  state 


_  P 

pa  =  Ra  «a 

where  is  the  gas  constant  for  dry  air,  and 
ture  in  degrees  Kelvin.  The  magnitude  of  the 
computed  at  each  grid  point  from 


(54) 

the  surface  air  tempera- 
geos  trophic  wind  was 


2.  Correction  for  Curvature 

A  gradient  wind  was  calculated  to  correct  for  the  effect  of 
curvature  on  the  surface  geostrophic  wind.  Two  basic  principles  were 
required  for  this  correction.  The  first  was  that  following  an  air  parcel 
along  an  isobar,  the  total  derivative  was  zero  or 


dP  (x , y)  =  0. 


(56) 


The  second  requirement  was  the  mathematical  statement  of  curvature 

Kc(x'y)  =  [i  +  (y)£"  ]  3/2  =  _R  (57) 

where  y’  and  y  *  re  the  first  and  second  derivatives  with  respect  to  x, 
and  R  was  the  radius  of  curvature .  By  expanding  the  total  derivative  of 
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P,  obtaining  y'  and  y",  and  rearranging,  the  following  result  was  ob- 

I 

talned  for  the  radius  of  curvature 


(_az,2 

n  &V 

,  ,  Sp/Sx  2 

L1+  J 

v2 

R  = 

&  P 
d  x 

a2  p  a2  p  ap/ax* 

ax  ay  1  a  y z  a  p/ay 

. 2  *  (58) 

ap  azp 

ay  ax2 

Next  the  gradient  wind  was  computed  from 


VG  =  G- 


V2 

G 

fR 


(59) 


Solving  for  V_  and  rearranging  led  to  a  workable  form  of  the  equation 
G 


= 


4- 


(60) 


G  !/2  +  U  1/4  +  G/Rf 

*  1 

For  this  case  the  magnitude  of  the  term  G/Rf  determined  how  much  curva- 

I 

ture  would  be  allowed.  The  maximum  values  were  chosen  to  be  -3/16 
for  anticyclonlc  curvature  and  +3/4  for  cyclonic  curvature  which  yielded 
a  range  of  =  1.33G  for  anticyclonic  curvature  and  =  0.67G  for 
cyclonic  curvature.  This  agreed  with  generally  accepted  synoptic  limits. 

Subroutine  GRADWD  (Appendix  A)  makes  the  necessary  com¬ 
putations  and  follows  a  similar  scheme  utilized  by  FNWC  (Kaitala,  1971 
unpublished).  Centered  finite  difference  schemes  were  used  to  calculate 
the  derivatives .  \ 

3 .  Thermal  Wind  Correction 

The  nondimensional  thermal  wind  components  were  deter¬ 


mined  fron1^ 


d  9 


(61) 


VT  = 


*2% 


b  9 

ax  ' 


(eh 


and  the  total  thermal  Wind  was  given  by 

r  =  y  u|  +  v|  . 


(63) 


For  this  model,  r  had  a  limiting  value  of  100.  The  geostrophic  and 
-.Iherraal  wind ^angles  were  calculated  from  the  previously  commuted  wirid 

components.  This  procedure  was  repeated  for  each  grid  point  and  was 

i  \ 

'  1  - 

used  as  input  data  to  the  boundary  layer  model. 

,  .  ■  ’  •  ] 

B.  PkOCEDURES  IN  SOLVING  THE  BOUNDARY  LAYER  EQUATIONS 

| 

1.  Initial  Conditions 

a.  \  Geostrophic  Wind  Speed  Less  Than  Or  Equal  to  5.0 
meters/second: 

This  case  was  associated  with  the  aerodynamically 

Smooth  flow  region  described  ^arlier.  For  this  case  n<!?  corrections  were 

]  \  ]  1 

made  for  stability  and  a  constant  drag  coefficient  of  0.022  was  assumed, 

This  resulted  in  a  constant  inflow  angle  of  15°.  The  19.5  meter  wind 


wa)s  determined  from 


\ 


U19.5  "  °‘7G 


(64) 


The  choice  of  the  drag  coefficient  w^s  not  arbitrary,  but  depended  on  the 

|  : 

relationship  between  C  and  U  (see  Figure  2).  The  value  of  Cin  was 

Z  Z  XU 

a  minimum  at  6  meters/secpnd  under  neutral  conditions.  The  value  of 

-  I  I 

Cjg  5  was  extrapolated  from  this  curve. 
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b.  Geostrophic  Wind  Speed  Greater  Than  5.0  meters/ 

second  and  Absolute  Value  of  Air-Sea  Temperature 

Difference  Less  Than  or  Equal  to  1°C: 

Initially  the  drag  coefficient  was  C  =  0.0245,  with 

z 

the  thermal  wind  (r)  and  the  angle  between  the  geostrophic  and  thermal 
wind  (  f| )  set  equal  to  zero.  For  this  case  L*,  L'  and  IJJ  (h/L’)  were 
set  equal  to  zero  and  0U  (h/L')=  1.0.  A  solution  was  then  obtained  using 
the  modified  barocllnic  model  of  Blackadar  described  below. 

c.  Geostrophic  Wind  Speed  Greater  Than  5.0  meters/ 
second  and  Absolute  Value  of  Air-Sea  Temperature 
Difference  Greater  than  1°C: 

The  drag  coefficient  C  =  0.0245  is  used  to  estimate 

z 

the  initial  U*  and  r  and  n  are  set  equal  to  zero.  &  -  0  was  restricted 

d  S 

to  the  range  -15  to  +4.0°C.  The  values  of  L#,  L',  III  (h/L'),  and 
jef  (h/L')  were  computed  from  z  ,  0  -0  and  U*  and  a  solution  was 

U  d  d  S 

obtained  using  the  modified  Blackadar  model  described  below. 

L#  was  computed  from  (34)  and  L'  from  (15)  using  as 

a  first  approximation 

UJ(za/L')  =  10  -  In  (za/zQ)  .  (65) 

III  (h/L')  was  computed  from  SUBROUTINE  PSI  (Appendix  A).  If  L#  was 
greater  than  zero  (stable  case),  then 

Ul(h/L')  =  0.7  L*  (66) 
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If  was  less  than  or  equal  to  zero  (unstable  case),  a  value  of 
Richardson  number  and  0U  (h/L')  had  to  be  computed.  This  was  accom¬ 
plished  in  SUBROUTINE  SHR  (Appendix  A).  Then  l|l  (h/L')  was  calculated 
from  (13). 

0u(h/L')  was  computed  using  SUBROUTINE  SHR 
(Panofsky,  1968;  Card  one,  1969).  If  L#  was  greater  than  zero  (stable 
case),  then 


0  (h/L1)  =  1  +  7L*  (67) 

u 

If  L*  was  less  than  or  equal  to  zero  (unstable  case),  then  a  new 
Richardson  number  was  determined  from 


R*new  ”  L* 
where  initially  Rj^  =  L+  . 
inequality 


1/4 

(1-18Rlold)  <68) 

This  iteration  was  continued  until  the 


Ri 

Jnew 

was  satisfied,  where 
found  from 


<  < 


(69) 


«  =  0.01  for  this  model.  Then  0U  (h/L1)  was 


-  ITTIeRl — ri/4  •  (7W 

*new  ' 

2.  The  Modified  Blackadar  Barocllnlc  Boundary  Laver  Model 
The  basis  of  this  model  was  an  Iteration  procedure  which 
Included  the  effects  of  stability  and  barocllniclty  on  the  computed 
wind.  The  procedure  was  free  of  computational  instability  and  converged 
rapidly  for  all  reasonable  choices  of  initial  U*  .  First  values  of  zQ  , 
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were  obtained  using  re ^ pec- 


P,  P\  »,  q,  *  ,  t  ,  UJ0,  and  U* 

°  new 

tlvely  equations  (23),  (35),  (38),  (43),  (44),  (45),  (46),  (50),  and  (51). 
The  Iterative  procedure  was  repeated  until  the  Inequality 

|U*  -  U*  I  S  <  (71) 

I  new 

was  satisfied  where  <  was  0.05.  The  Iteration  was  repeated,  setting 

U*«U*  ,  until  the  Inequality  was  again  satisfied.  If  L*  was  not 

new 

zero,  then  new  values  of  Ul  ,  zQ  ,  and  L'new  were  computed  where  the 
argument  of  Ul  was  za/L'.  If 

L'„ew-L'  *  5  <”> 

then  L*  was  set  equal  to  L'new  .  The  model  was  repeated  from  where 
L*  was  first  calculated  until  the  inequality  was  satisfied.  After  com¬ 
pleting  the  Iterations  for  the  case  where  the  thermal  wind  (r)  and  I) 
were  equal  to  zero,  the  computed  values  of  r  and  n  were  substituted 
and  the  model  was  again  repeated  until  all  criteria  were  satisfied.  If 
LVO,  then0(z/L’)  was  computed  as  0(19-  5/L* )  ;  otherwise  0(19. 5/L* ) 
was  set  equal  to  zero.  Finally  the  wind  at  19.5  meters  was  determined 
from 


Values  of  0fl - .  G,  zq  ,  z/L1  ,  ,  U#  ,  UJ9  $  ,  log  Rq  ,  U^/G, 

Ujg  g/G  and  Rj  were  then  available  for  each  grid  point. 
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C.  COMPUTATION  OF  INDEX  OF  REFRACTION  CONSTANT 


In  this  section  possible  relationships  are  considered  between  the 
refractive  index  structure  constant  and  quantities  dependent  only  on  the 
wind  speed  and  air-sea  temperature  difference.  Specifically  those 
quantities  that  can  be  predicted  from  the  model  considered  In  this  study 
are  examined. 

The  propagation  of  light  waves  and  radio  waves  Is  affected  pri¬ 
marily  by  air-density  inhomogeneltles  and  therefore  by  temperature  fluc¬ 
tuations.  One  effect  of  these  fluctuations  is  scintillation  or  rapid  changes 
in  intensity  at  a  point.  These  are  distortions  caused  by  variations  In  the 
speed  of  propagation  through  turbulent  clemonts  having  different  densi¬ 
ties,  and  hence  different  Indexes  of  refraction.  The  rofractive  Index 
structure  constant  (C^  Is  the  primary  parameter  characterizing  the  turbu¬ 
lent  f’uctuatlons  of  the  atmospheric  t  jfractlve  Index  for  scales  from  a 

2 

few  millimeters  to  several  meters.  If  the  quantity  Cn  Is  known,  cal¬ 
culations  can  be  made  on  scintillation,  beam  spreading,  and  beam  wander, 
all  of  which  arc  Important  in  describing  light  wave  propagation. 

A  relation  necessary  in  the  formulation  is  the  one-dimensional 
spectrum  for  temperature  fluctuations  in  the  inertial  subrange  described 
by  CorTsln  (1951) 

-  V3  —Vi 

0g  00  ■  B  C  N  K  (75) 

where  c  is  the  rate  of  dissipation  of  turbulent  kinetic  energy,  N  is  a 
moasure  of  the  “smearing"  of  temperature  inhomogeneities.  B  is  an 
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enpirlcal  constant  ("0.01,  Boston,  1970),  and  Kli  tha  streamwise 

component  of  wave  number.  The  structure  constant  for  temperature  flue- 
2 

tuatlons  (Cj)  Is  defined  as 

2  *1/ 

CT  -  B  C  N  (76) 

and  Is  related  to  the  Index  of  refraction  as  follows 

C2 .  d  .  (77) 

n  ^  T 

A 

Another  expression  which  has  to  be  considered  Is  tho  nondlmon- 
slonal  balance  equation  for  turbulent  kinetic  energy,  which,  for  station¬ 
ary,  horizontally  homogeneous  flow  Is  (Busch  and  Panofsky,  1968) 


-  */l  -  K  -  *T 


o. 


(78) 


0U  Is  the  nondlmenslonal  wind  shear  (8),  0D  represents  toe  diver* 
gence  terms ,  and 


ikt 

U? 


(79) 


The  following  rotations  arise  from  the  definitions  in  tho  non- 


dimensional  turbulont  kinetic  energy  balance  and  definition  of  N 


kz 


(80) 


N  ■  -«w  - 


(81) 


Substitution  of  these  two  expressions,  plus  the  definition  of  K  (29) 

m 


into  the  definition  of  C*  yields  (Panofsky,  1968) 

<4  * B *t/3  'ii’2- 


(82) 
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Equation  (82)  la  not  a  aultabla  expression  because  of  the  lack  of 
Information  on  0t  However  If  the  effect  of  the  divergence  terms 

Is  neglected,  can  be  expressed  in  terms  of  0 and  t/L  as 


0t  •  0U  -  f/L  .  (83) 

It  Is  noted  that  with  wind-wave  Interaction,  the  assumption  that  the 

divergence  terms  can  be  neglected  Is  probably  not  valid,  but  the  possible 

2 

effect  this  assumption  has  cn  predictions  of  C^.  has  not  been  studlod. 


Substitution  of  (83)  Into  (82)  ylolds 


r1/ 3 


-  B  0  - t/L 
u 


1 

Km  *u 
»  ■ 


4/3  2 
(kz)  C-ff")-  (84) 


This  final  equation  for  C2  lends  Itself  to  the  boundary  layer  model  be- 

T 

cause  all  the  parameters  can  be  computed  within  the  model.  Thus,  the 

2 

model  can  be  used  to  obtain  a  horizontal  distribution  of  over  water 
for  various  wind  and  stability  conditions.  The  average  stability  d9/dz 
Is  determined  from  (14). 

The  ratio  of  the  exchange  coefficients  K^/K^  U  sot  equal  to  1.0 
In  Cardone's  formulation.  However  It  has  been  shown  that  and 
are  not  usually  Identical  except  for  small  values  of  d9/  ds  for  which 
the  turbulent  motion  Is  not  significantly  influenced  by  buoyancy  forcos. 
Businger  et  al.  (1971)  suggested  the  following  form  for  the  ratio  of  eddy 


diffusions 


•h  -  1.35 


(1-9  VL)  /2 
( 1  -  1 S  */L)  */4 
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Hero  a.  *  1.35  for  neutral  condition!,  which  agrees  closely  with  labora- 
n 

tory  maaturamants  (c(. ,  Hinxe,  19S9).  A  valua  of  1.3S  rathar  than  1.0 
for  a^  undar  nautral  conditions  arliai  bacausa  tha  von  KArroAn  constant 
was  equal  to  0.3S  Instead  of  tha  usual  0.40  (Buslnger  at  at. ,  1971). 


D.  WAVE  HEIGHT  ANALYSIS 

Wava  analysis  Is  another  application  for  tho  boundary  layer  modal. 
Tor  this  purpose,  a  simplified  version  of  tha  TNWC  sea-swell  model 
(Hubert  and  Mendenhall,  1970)  was  utilised.  The  u  and  v  components 
of  the  19. 5  meter  wind  determined  from  the  boundary  layer  model  were 
used  as  both  current  and  history  Input  In  this  analysis. 

Tha  TNWC  sea-swell  modol  computes  a  weighted  wind  from  tha 
wind  components  as  follows: 

u  -  0.3ut_J2  ♦  0.7  ut 


v  •  0.3vt_j2  ♦  0.7Vj 


(86) 


where  t- 12  Is  the  12-hour  old  history  wind  speed  and  t  Is  the  current 
wind  speed.  Those  component  values  are  saved  as  history  for  the  next 
analysis  run.  and  then  the  total  wind  (V)  is  computed.  The  wave  heights 
(WA )  in  feet  are  then  computed  from 


7 


(87) 


Although  this  Is  a  very  simple  empirical  method  for  determining  wave 
heights,  the  basic  method  has  been  used  by  TNWC  since  the  mid- 1960* s 
and  this  particular  model  has  been  used  operationally  since  August  1970. 
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In  addition,  tha  FNWC  modal  computad  corractlona  for  fatch  and 
lea,  and  parformad  soma  smoothing  of  tha  raaultlng  halght  flald.  It  was 
assumad  thasa  corrections  would  ba  minimal  for  this  application. 
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.  IV.  SPECIFICATION  OF  EXTERNAL  PARAMETERS 


A  primary  objective  of  this  study  was  to  apply  the  two-level  baro- 
cllnlc  boundary  layer  model  to  actual  atmospheric  conditions  and  compare 
the  computed  results  with  the  observed  results.  The  time  period  00  GMT 
28  November  69  through  00  GMT  02  December  69,  with  an  analysis  Inter¬ 
val  of  every  12  hours,  was  chosen  for  this  comparison  for  several  reasons. 
First  an  intense  low  pressure  system,  which  influenced  most  of  the  central 
North  Pacific  Ocean,  developed  during  this  period  with  the  system  reach¬ 
ing  maximum  intensity  about  00  GMT  30  November  69 .  The  system  pro¬ 
vided  Input  data  which  allowed  the  model  to  predict  boundary  layer 
parameters  under  varying  atmospheric  conditions.  Secondly  considerable 
synoptic  data  had  been  gathered  by  the  author  during  an  earlier  study  on 
the  extreme  sea  conditions  (unpublished)  associated  with  this  storm.  The 
latter  gave  a  good  basis  for  comparing  computed  values  to  observed  data. 
Finally,  this  storm  produced  waves  and  swell  which  caused  considerable 
damage  to  several  Islands  of  the  Pacific  as  well  as  to  the  spar  buoy  FLIP. 
Therefore  a  wave  generation  model  was  applied  using  the  predicted  sur¬ 
face  winds,  and  these  predicted  waves  were  compared  to  the  observed 
wave  conditions . 

The  00  GMT  30  November  69  period  was  of  the  most  interest  due 
to  Its  extreme  conditions,  extensive  data  coverage,  and  it  was  the  prob¬ 
able  period  when  the  destructive  waves  were  being  generated.  The  various 
corrections  described  above  were  incorporated  into  the  model  for  this 
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time  period ,  and  synoptic  data  were  plotted  and  analyzed  for  surface 
pressure,  air-sea  temperature  difference,  surface  wind  speed  and  wave 
height. 

Surface  layer  parameters  were  computed  Including  various  comblna 
tions  of  stability,  curvature,  and  thermal  wind  to  yield  Information  on 
the  relative  Importance  of  each  parameter  In  the  formulation.  The  follow 
ing  notation  will  be  used  In  describing  the  results: 

NN  -  non-neutral,  correction  for  stability 
N  -  neutral  -  no  stability  correction 
NC  -  no  curvature  correction  applied 
C  -  curvature  correction  applied 
NT  -  no  thermal  wind  correction  applied 
T  -  thermal  wind  correction  applied 

It  was  recognized  at  the  onset  of  the  study  that  the  specification 
of  the  external  fields  (  A  9  and  G)  had  to  be  reasonably  accurate  In 
order  to  evaluate  the  model  on  the  basis  of  observed  wind  fields.  For 
example,  if  the  observed  surface  winds  were  quite  different  than  those 
predicted  by  the  model,  some  estimate  had  to  be  available  on  the  possi¬ 
bility  that  the  discrepancies  were  due  to  poor  specification  of  the  geo- 
stophlc  wind  or  the  air-sea  temperature  difference.  Therefore  In  this 
section  consideration  is  given  to  the  relevant  aspects  of  representing 
those  fields  which  make  up  the  external  parameters . 

The  charts  in  the  following  figures  correspond  to  a  1:30  million 
polar  stereographic  projection  with  the  grid  as  depicted  in  Figure  9 . 
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A.  FNWC  ANALYSIS 


Sea  level  pressure  and  air-sea  temperature  difference  fields  for 
00  GMT  30  November  appear  in  Figure  10.  A  967  mb  low  pressure 
center  was  located  near  45 °N  180°  with  considerable  thermal  instability 
both  ahead  of  and  behind  the  system.  A  significant  feature  of  this  sys¬ 
tem  was  the  intense  isobaric  gradient  in  the  western  quadrant. 

The  first  consideration  was  the  comparison  of  the  computed  surface 
geostrophlc  wind  with  the  computed  gradient  wind  for  this  period.  The 
geostrophic  and  gradient  wind  fields  computed  from  the  surface  pressure 
fields  in  Figure  10  appear  in  Figures  11  and  12  (units  of  meters/second). 
Figure  13  depicts  the  Influence  of  the  curvature  correction  on  the  geo¬ 
strophic  wind  computed  as  a  percent  difference  from  geostrophic  wind 
minus  gradient  wind  then  divided  by  geostrophic  wind . 

The  two  maximum  wind  centers  west  and  north  of  the  Low  are  based 
on  data  at  one  grid  point  and  result  from  the  influence  of  anticyclonic 
curvature,  strong  isobaric  gradients,  and  finite  difference  computations. 
These  features  do  not  appear  to  be  realistic.  However,  the  wind  speed 
distribution  depicted  by  them  is  good  except  that  the  magnitude  of  the 
wind  speed  may  be  a  little  too  large.  It  was  expected  that  the  maximum 
change  would  occur  in  an  area  of  maximum  curvature  and  maximum  pres¬ 
sure  gradient  from  the  relation  in  (60).  Generally  speaking,  the  change 
in  wind  speed  due  to  the  curvature  correction  is  10  to  20%  with  a  maxi¬ 
mum  change  of  34%  imposed  by  the  limits  placed  on  the  term  G/Rf  in  (60). 
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Geoatrophlc  wind  speed  analysis  based  on  surface 
pressure  (Figure  10) 
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Gradient  wind  speed  analysis  based  on  surfaoe 
pressure  (Figure  10) 


B.  SUBJECTIVE  ANALYSIS 


The  results  from  the  hand  analyzed  surface  chart  were  also  used 
In  comparisons  of  results.  Figure  14  Illustrates  the  surface  pressure 
and  air-sea  temperature  difference  fields  from  the  subjective  analysis 
which  are  analogous  to  Figure  10.  The  major  variations  between  Figure 
10  and  14  appear  to  be  due  to  the  difference  In  the  number  of  synoptic 
reports  used  as  well  as  the  analysis  scheme.  Theoretically,  Figure  14 
should  be  more  representative  since  more  data  were  included,  but  it  was 
a  subjective  analysis.  Figures  IS  and  16  depict  the  geostrophic  and 
gradient  surface  winds  computed  from  the  subjectively  analyzed  pressure 


field. 
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Surface  pressure  (solid)  and  sir- tea  ten pars 
(dotted)  fields  for  00  OCT  30  Mot  69  based  a 
plotted  data 
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Or&dlont  wind  >r«*d  flold  taood  on  aurfaeo 
rrtuurc  (Firuro  14) 


V.  RESULTS  FROM  EXPERIMENTS  WITH  THE  MODEL 


A.  COMPARISONS  IN  TERMS  OF  SURFACE  ROSSBY  NUMBER  AND 
STABILITY  CATEGORIES 

Near  surface  parameters  were  computed  for  the  19.5  meter  level. 

This  level  represents  an  average  height  of  ship  anemometers  (Bunting, 

1968).  Theoretically  the  computed  19.5  meter  level  winds  could  be 

compared  with  the  observed  wind  fields  and  this  would  yield  some 

information  on  the  accuracy  of  the  model  in  predicting  surface  winds . 

In  order  to  verify  the  computations  and  also  to  examine  the  range 

of  conditions  provided  by  the  data,  a  comparison  was  made  with  Cardone's 

results  which  appear  in  Figures  4  and  5.  Results  for  this  comparison 

were  obtained  from  computations  in  which  stability  was  included  but  not 

baroclinicity.  Results  from  all  nine  time  periods  appear  in  Figures  17 

and  18  where  values  of  the  logarithm  of  the  surface  Rossby  number  (24) 

range  from  7.76  to  9.45.  The  U./G  versus  log  R  results  (Figure  17) 

*  o 

agree  with  Cardone’s  results  (Figure  4)  except  with  regards  to  the  inter¬ 
val  between  L*  isolines.  In  Figure  17  this  interval  (spacing)  is  larger. 
With  respect  to  versus  log  RQ  (Figure  18),  the  unstable  results 
(negative  L*)  agree  with  those  in  Figure  5  except,  again,  with  regards 
to  the  interval  between  L*  isolines.  The  stable  results  (positive  L*), 
however,  deviate  considerably  from  Cardone's  results  with  the  veering 
angle  decreasing  as  L*  increases.  The  results  in  Figure  18  would 
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indicate  that  the  veering  angle  decreaiei  with  both  Increasing  stability 
apd  increasing  Instability  with  a  discontinuity  occurring  near  neutral 
conditions . 

The  unusual,  and  unexpected,  features  in  the  results  Just  de¬ 
scribed  were  examined  In  considerable  detail.  The  same  results  were 
noted  by  Paflas  (unpublished  study,  1971)  in  an  Investigation  with 
Cardone's  program.  These  features  in  the  results  are  perhaps  related 
to  the  drag  coefficient  representation  selected  by  Cardone  and  also  be¬ 
cause  different  iteration  limits  were  used  in  this  study  than  those  used 
by  Cardone. 

The  range  of  data  with  respect  to  U^/G  and  0^-0^  was  also 
examined  and  the  results  appear  in  Figure  19.  The  computations  for 
these  results  were  made  on  the  basis  of  a  barotropic  atmosphere  with  no 
correction  for  lsobaric  curvature.  These  computations  differ  from  those 
used  by  Cardone  for  Figure  8  in  that  a  variable  coriolis  parameter  was 
used.  Isolines  have  been  drawn  for  geostrophic  wind  speeds  ranging 
from  15  to  75  knots.  For  those  wind  speeds  considered  by  Cardone  in 
Figure  19,  and  also  for  wind  speeds  above  20  knots  considered  by 
Paflas,  these  results  are  in  agreement.  The  change  in  shape  of  the 
curves  for  winds  under  20  knots  was  not  described  by  Cardone  but  is 
probably  due  to  a  combination  of  the  functional  forms  for  the  stability 
dependence  and  the  drag  coefficient  as  chosen  by  him. 
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Fig.  ^9.  Geostrophlc  drag  coefficient  versus  alr~sea 

temperature  difference  for  range  of  geostrophlc  wind 
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B.  EXPERIMENTS  WITH  CHANGING  THE  PHYSICS  IN  THE  TWO- LAYER 

MODEL  AND  CHANGING  THE  SPECIFICATION  OF  THE  SURFACE 

GEOSTROPHIC  WIND 

In  this  section  U^g  results  which  were  obtained  from  different 
computations  with  the  two-layer  model  are  compared.  These  computations 
were  varied  by  including  or  not  including  specific  parameters  or  correc¬ 
tions  .  The  variations  considered  and  compared  were  described  earlier 
and  are  1)  N/NC/NT,  2)  N N/NC/NT,  3)  NN/NCA,  4)  NN/CA,  and 
5)  NN/NCA  with  k=  0.35. 

These  variations  will  be  referred  to  as  Cases  1  through  5.  It  is 
noted  that  if  the  geostrophic  wind  speed  was  less  than  or  equal  to  5.0 
meters/second ,  stability  and  baroclinicity  were  not  taken  into  account 
in  the  iteration  schemes .  In  regions  with  this  wind  speed  range ,  espe¬ 
cially  in  the  subtropics,  differences  will  not  arise  except  where  isobarlc 
curvature  exists. 

Per  cent  difference  (1-2)  will  be  used  to  compare  the  two  cases 
where  (1-2)  corresponds  to  the  following  computation  at  each  grid  point: 

case  1  result  -  case  2  result  Y  1fjn 
case  1  result 

Although  large  differences  between  cases  were  unexpected,  small  cor¬ 
rections  at  low  wind  speeds  could  produce  large  per  cent  changes. 

The  surface  wind  fields  corresponding  to  each  of  the  five  cases 
are  presented  in  the  following  paragraphs  along  with  selected  comparisons 
between  pairs  of  cases . 
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Case  1  (N/NC/NT) 

This  Is  the  simplest  case  with  respect  to  the  possible  physics 
which  are  Included  In  the  model.  The  external  wind  field  is  defined  for 
a  neutral,  barotropic  atmosphere  with  no  curvature  correction.  The  wind 
field  at  19.5  meters  for  this  case  appears  in  Figure  20. 

...  „  .  Case  2  (NN/NC/NT) 

The  effects  of  stability  are  now  included  and  the  19 . 5  meter  wind 
field  for  this  case  appears  in  Figure  21.  The  per  cent  difference  (2-1) 
field  is  shown  in  Figure  22.  Positive  values  correspond  to  increases  in 
the  computed  U1Q  _  field  due  to  including  the  stability  effect.  This  in- 
crease  is  quite  evident  on  the  western  side  of  the  cyclone.  The  patterns 
are  not  expected  to  be  coincident  with  the  air-water  temperature  isolines 
in  Figure  10  because  the  stability  effect  is  dependent  on  both  the  temper¬ 
ature  gradient  and  the  wind  speed.  The  20%  contours  near  the  eastern 
border  are  attributed  to  small  corrections  at  low  wind  speeds  as  dis¬ 
cussed  previously. 

Case  3  (NN/NC/T) 

The  effects  of  including  the  baroclinicity  or  thermal  wind  are  now 
included  and  the  19.5  meter  wind  field  for  this  case  appears  in  Figure 
23.  The  per  cent  difference  (2-3)  field  is  seen  in  Figure  24.  Negative 
values  correspond  to  an  increase  in  the  computed  19.5  meter  wind  due 
to  including  the  thermal  wind  in  the  model.  In  this  difference  field,  only 
20  grid  points  (6%)  appeared  to  be  significantly  affected  by  including  the 
thermal  wind .  These  grid  points  generally  lie  in  an  area  of  warm  air 


advectlon. 
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19*5  meter  wind  difference  field  (%)  for  Case  2  < 
Case  3.  Negative  values  indicate  increased  wind 
arjoed  due  to  thermal  wind 


£>11 4  (W/g/T) 

In  addition  to  all  specifications  possible  within  the  physics  of  the 
model,  tho  correction  for  curvature  Is  applied  to  the  external  wind  flold. 
The  19. 5  meter  wind  analysis  for  this  field  Is  shown  in  Figure  25  and 
the  per  cent  difference  field  (3-4)  appears  in  Figure  26.  Negative  values 
correspond  to  larger  19.5  meter  winds  due  to  including  the  curvature  cor¬ 
rection.  Strong  wind  maxlmums,  noted  in  comparing  the  external  wind 
fields  with  and  without  this  correction,  also  appear  in  Figure  26  as  the 
predominant  feature.  In  general,  the  effective  change  in  the  19.5  meter 
wind  field  specification  is  larger  and  covers  more  area  t'lan  any  of  the 
other  cases . 

Case  5  (NN/NC/T  with  k«0.35) 

This  case  differs  from  case  3  only  with  respect  to  the  change  of 
the  von  Karmen  constant  from  0.40  to  0.35.  The  19.5  meter  wind  field 
for  this  case  appears  in  Figure  27  and  the  per  cent  difference  field  (3-5) 
is  depicted  in  Figure  28.  Negative  values  correspond  to  increased  19.5 
meter  winds  due  to  the  change  In  the  value  of  k.  The  suggested  value 
of  k=0.35  arose  from  recent  surface  layer  observational  investigations 
by  Buslnger  et  al.  (1971).  They  noted  from  their  results  that  using 
k=0.40  in  several  empirical  expressions  led  to  estimates  of  U*  values 
which  were  15%  too  high  when  compared  with  direct  measurements  of 
the  momentum  transfer.  The  fields  in  Figure  28  indicate  that  a  general 
increase  of  12  to  17%  occurs  in  the  computed  19.5  meter  wind  field  if 
the  lower  value  of  k  is  used.  This  increase  agrees  with  the  U* 
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meter  wind  speed  analysis  (m. 
4  (NN/C/T)  based  on  FNtfC  data 
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discrepancies  noted  by  Businger  et  al. .  Because  the  von  Klrmtfn 
constant  appears  In  several  expressions  In  the  boundary  layer  model, 
the  Importance  of  its  change  In  the  two-layer  model  was  not  known. 
These  results  indicate  that  the  difference  between  k- 0.35  and  k*  0.40 
is  important. 

In  summary,  the  five  cases  described  above  indicate  which  terms 
were  most  influential  in  the  model.  It  should  be  remembered  that  most 
of  the  subtropical  point?  were  not  Influenced  by  the  model  since  the 
wind  speeds  were  generally  less  than  5.0  meters/second.  Of  all  the 
modifications  to  the  model,  the  curvature  correction  to  the  external 
field  produced  the  largest  variation  over  the  entire  grid .  Changing  k 
from  0.40  to  0.35  accounted  for  the  second  most  significant  change. 
Differences  in  the  surface  wind  fields  were  evident,  however,  when 
specifications  in  the  model  were  made.  Of  the  two  possible  specifica¬ 
tions,  stability  appeared  to  be  more  significant  than  the  thermal  wind 
with  respect  to  area  influenced  and  per  cent  difference  changes. 

Although  the  influence  of  the  thermal  wind  in  this  model  appears 
to  be  small  with  respect  to  the  relative  magnitude  of  the  computed  wind, 
the  influence  it  has  on  the  direction  of  the  wind  appears  to  be  signifi¬ 
cant.  The  difference  field,  in  degrees,  of  the  wind  direction  for 
NN/NCA  -  NN/NC/NT  appears  in  Figure  29.  This  field  represents  the 
difference  in  wind  direction  between  a  case  with  and  a  case  without  the 
thermal  wind  specification.  The  region  of  significant  differences  does 
not  coincide  with  that  region  where  the  thermal  wind  effect  was  large 
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19.5  meter  wind  direction  difference  fleld(deg 
for  Caae  3  -  Case  2.  Positive  values  Indicate 
Increased  or  veering  wind  direction 


with  respect  to  wind  speed.  Positive  values  in  Figure  29  correspond  to 
a  veering  wind  and  negative  values  correspond  to  a  backing  wind.  It  is 
noted  that  the  direction  change  due  to  baroclinicity  is  primarily  deter¬ 
mined  by  the  angle  between  the  thermal  and  geostrophic  winds . 

C.  COMPARISONS  OF  COMPUTED  WINDS  WITH  OBSERVED  WINDS 

These  comparisons  will  be  described  for  the  case  NN/C/T  with 
k=0.35.  The  19.5  meter  wind  fields  computed  from  the  other  cases  were 
compared  to  the  subjectively  analyzed  observed  surface  wind  field 
(Figure  30)  and  in  all  cases  the  computed  winds  underestimate  the  ob¬ 
served  winds.  This  was  also  the  case  when  the  above  corrections  were 
included,  but  yielded  the  best  results  with  respect  to  the  magnitude  and 
patterns  of  the  wind  speeds . 

Because  the  previous  comparisons  revealed  the  importance  of  the 
specifications  of  the  external  wind  field  with  respect  to  curvature,  it 
was  considered  necessary  in  this  comparison  to  examine  the  importance 
of  the  Initial  analysis  (FNWC  or  subjective).  The  19.5  meter  wind  field 
computed  with  the  above  specifications  and  with  the  FNWC  analysis  is 
shown  in  Figure  31  and  the  wind  field  based  on  subjective  analysis  is 
seen  in  Figure  32.  The  observed  surface  wind  field  which  appears  in 
Figure  30  was  analyzed  using  only  those  observations  reported  above 
10  meters/second.  In  Figures  31  and  32,  the  maximum  contour,  30 
meters/second,  agrees  with  the  observed  maximum  wind  speed,  but  the 
objective  analysis  does  not  indicate  the  25  meter/second  maximum  east 

of  the  low . 
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A  comparison  of  computed  19 . 5  meter  winds  from  the  two  analyses 
(Figures  31  and  32)  and  the  external  wind  fields  (Figures  11  and  15) 
shows  that  differences  in  the  external  wind  fields  are  consistently 
carried  through  the  model.  The  model,  however,  was  able  to  produce 
19.5  meter  wind  fields  in  both  cases  which  were  in  good  agreement  with 
the  observed  wind  field,  at  least  for  wind  speeds  over  10  meters/second. 

D.  RESULTS  FROM  WAVE  HEIGHT  COMPUTATIONS 

Wave  height  fields  were  computed  using  the  19.5  meter  wind 
fields  obtained  from  the  boundary  layer  model  and  the  simplified  wave 
height  program  developed  from  the  FNWC  operational  sea-swell  model. 
The  program  was  run  starting  at  00  GMT  28  Nov  in  order  to  build  history 
for  this  analysis.  The  19.5  meter  wind  computations  were  obtained 
using  all  specifications  and  corrections  described  previously,  and  they  . 
were  also  made  on  both  the  FNWC  and  subjectively  analyzed  fields.  The 
computed  wave  field  from  the  FNWC  data  is  seen  in  Figure  33  and  the 
wave  field  from  the  subjective  analysis  appears  in  Figure  34.  The  ob¬ 
served  wave  field  obtained  from  synoptic  observations  is  depicted  in 
Figure  35.  It  is  noted  that  wave  height  is  probably  the  most  difficult 
parameter  to  observe  and  therefore  a  representative  analysis  of  the  ob¬ 
served  wave  field  is  difficult  to  achieve. 

The  wave  field  obtained  from  the  subjective  analysis,  Figure  34, 
appears  to  be  in  good  agreement  with  the  observed  heights  with  maxi¬ 
mum  computed  heights  slightly  less  than  the  maximum  observed  heights. 
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Computed  wave  height  analysis  (feet)  based  on 
data  for  00  GMT  30  Nov  69 
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Observed  surface  wave  height  analysis  (feet)  based 
on  plotted  data  and  subjectively  analyzed 


Results  computed  from  the  FNWC  analysis,  Figure  33,  appear  to  have 
similar  agreement  with  respect  to  the  maximum  height  but  poorer  agree¬ 
ment  with  respect  to  the  overall  pattern  of  the  wave  field. 

In  general  these  results  indicate  that  the  two-layer  model  can 
produce  surface  wind  fields  which  can  be  used  for  wave  height  analysis, 
even  from  a  very  simple  wave  model.  It  is  noted  that  in  the  subtropics 
the  wave  heights  computed  from  both  sets  of  analysis  are  considerably 
below  the  observed  wave  heights.  This  is  probably  due  tolow  surface 
geostrophic  wind  speeds  in  this  area  and  the  fact  that  no  corrections  are 
made  for  stability  or  baroclinicity  at  low  wind  speeds.  This  could  pos¬ 
sibly  be  corrected  in  the  model  by  selecting  a  higher  initial  value  for 
the  drag  coefficient. 

E.  RESULTS  FROM  INDEX  OF  REFRACTION  COMPUTATIONS 

A  final  application  of  the  boundary  layer  model  was  the  computa¬ 
tion  of  the  index  of  refraction  structure  constant  (C^)  using  (84).  The 
ratio  of  the  exchange  coefficients  was  determined  from  (85)  and  the 
input  data  for  this  equation  were  obtained  by  including  the  stability 
effect,  the  thermal  wind  effect,  the  correction  for  curvature,  and  using 
k=0.35.  All  parameters  needed  to  compute  C.2  were  determined  in  the 
model. 

The  results  from  these  computations  appear  in  Figure  36.  Values 

2 

of  range  from  0.0  in  regions  of  near  neutral  stratification  to  a  maxi¬ 
mum  of  0.181  in  the  region  of  maximum  wind  speeds.  Comparisons  of 
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these  values  with  those  measured  by  Portman  et  al.  (1968)  Indicate  that 
they  are  comparable  in  regions  of  low  wind  speed  with  0-0  less  than 

d  S 

2 

1  degree.  Although  no  definite  relation  between  and  wind  speed  or 

2 

0-0  appears  in  Figure  36,  larger  values  of  C_  appear  to  be  associ- 
a  s  i 

ated  with  the  greater  wind  speeds  and  larger  instability  (0-0  ranging 

3  S 

.  .  from  -1  to  -6°C). 

2 

The  results  indicate  that  a  horizontal  distribution  of  C^,  for  a 
given  height,  in  this  case  19.5  meters,  can  be  computed  with  the  model 
with  some  degree  of  confidence  in  the  range  of  values  obtained.  Com¬ 
parisons  should  be  made  between  results  computed  by  this  model  and 
2 

observations  of  CT  to  obtain  better  verification. 
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VI.  CONCLUSIONS 


The  results  of  this  study  have  indicated  that  a  marine  boundary 
layer  model  can  be  applied  to  synoptic  scale  data  to  obtain  various  sur¬ 
face  parameters.  These  parameters  include  friction  velocity  (U*), 
surface  roughness  (zq),  surface  Rossby  number  (Rq),  19 . 5  meter  wind 
(Ujg  5),  veering  angle  (4>0)/  and  stability  length  (z/L). 

First,  the  results  Indicate  that  the  curvature  correction  produced 
the  greatest  change  to  the  wind  analysis.  However,  stability  and  modi¬ 
fication  of  the  von  Karman  constant  produced  significant  changes  in  the 
field,  while  changes  due  to  the  thermal  wind  were  confined  primarily  to 
changes  in  wind  direction. 

Second,  the  19.5  meter  wind  was  computed  including  various 
physical  specifications  and  corrections  and  compared  with  synoptic 
observations.  The  results  indicate  that  the  computed  wind  was  lower 
than  the  observed  wind  in  all  areas  and  for  all  cases,  but  the  general 
patterns  agreed  well  with  observations. 

Finally,  it  appears  that  Cardone's  planetary  boundary  model  pro¬ 
vides  a  reasonable  framework  to  specify  those  surface  layer  parameters 
required  for  spectral  wave  forecasting  and  index  of  refraction  calculations. 

Further  studies  are  needed  to  more  completely  verify  and  improve 
the  model.  The  specification  of  the  drag  coefficient  is  probably  the 
aspect  in  this  model  which  needs  the  most  attention.  The  drag  coefficient 
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has  received  considerable  examination  in  observational  Investigations 

and  the  results  indicate  some  disagreement  among  investigators .  A 

suggestion  of  this  study  is  to  increase  the  initial  value  of  C  to  increase 

z 

the  computed  surface  wind  speeds  especially  in  areas  of  low  wind  speed. 
Application  of  this  model  to  a  region  of  reliable  synoptic  data  such  as 
during  the  BOMEX  experiment  in  1969  would  reduce  some  of  the  errors 
due  to  data  input  and  also  test  the  physics  of  the  model  in  the  tropical 
atmosphere.  Finally,  observations  of  the  index  of  refraction  over  a  wide 
range  of  wind  speed  and  stability  would  allow  reliable  comparisons  to  be 
made  of  the  structure  constant  (C-j>). 


PROGRAM  TO  COMPUTE  A  WIND  AT  19.5  METERS  USING  BLACKADAR  BARCCLINIC 
MODEL  OVER  A  SPECIFIED  GRID 
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COMPUTE  GEQSTROPH I C  WIND  FROM  SFC  P  FIELD 
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ABSOLUTS  VALUE  OF  AIR-SEA  TEMPERATURE  DIFFERENCE  LE 
CONSTANT  DRAG  COEFFICIENT  =  0,0245 
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